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1 .  INTRODUCTION 


The  problem  of  numerical  grid  generation  is  of  current  interest  in  many 
branches  of  engineering  particularly  in  aeronautics,  mechanical,  and  civil 
engineering.  The  spatial  grids  are  generated  either  by  algebraic  methods 
using  various  spline  and  transfinite  interpolations  or  by  solving  certain 
partial  differential  equations.  In  this  regard,  a  book  [l],  review  articles 
[2],  [3],  and  conference  proceedings  [JJ,  5,  6]  may  be  consulted. 

The  research  accomplishments  noted  herein  are  in  reference  to  the 
problem  of  grid  generation  in  a  given  surface  by  solving  a  set  of  elliptic 
partial  differential  equations.  The  mathematical  model  used  in  the  present 
research  has  been  developed  by  the  author  under  this  and  previous  AFOSR 
grants  and  has  been  described  in  publications  [7  -  1 6 J .  Reference  [16]  has 
been  attached  as  an  appendix  in  which  Eqs.  (19)  and  (21)  describe  the  mathe¬ 
matical  model. 

The  developed  mathematical  model  has  been  used  to  generate  the  Car¬ 
tesian  coordinates  as  functions  of  the  curvilinear  coordinates  when  the 
surface  in  which  the  coordinates  are  to  be  generated  has  been  specified 
either  analytically  or  by  discrete  data  points.  In  most  practical  situa¬ 
tions  the  surface  is  usually  specified  by  discrete  data,  and  therefore, 
there  was  a  need  to  develop  computer  routines  to  fit  a  global  equation  of 
the  form  F(x,y,z)  =  0  to  describe  the  surface.  For  complicated  body  shapes 
e.g.  an  airplane,  a  segmentation  technique  has  been  used  in  which  the 
surface  is  divided  into  suitable  sections  and  then  for  each  segment  the 
function  is  separately  generated.  The  need  for  specifying  the  function  F  is 
due  to  the  fact  that  the  forcing  term  R  in  Eqs.  (19)  or  (21)  of  the  appendix 
is  obtained  from  F.  In  the  period  under  consideration  we  have  devised  vari- 
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ous  routines  for  the  above  noted  purpose  and  they  are  described  in  Section 
2.  The  developed  computer  code  has  been  used  in  many  complicated  shaped 
surfaces.  Figure  1  shows  some  current  results.  Other  results  are  available 
in  the  published  cited  references. 

2.  NUMERICAL  SCHEMES 

A  computer  program  for  the  numerical  solution  of  Eq.  (21  )  has  been 
developed  by  using  point  and  line  SOR .  The  essential  difference  between  the 
grid  generation  in  a  flat  space  and  in  a  curved  surface  is  in  the  appearance 
of  the  forcing  term  on  the  right  of  Eqs.  (19)  and  (21).  The  term  R  in  Eqs. 
(19)  and  (21)  depends  on  the  first  and  second  partial  derivatives  of  the 
fitted  function  F  with  respect  to  the  Cartesian  coordinates.  (Refer  to  [15] 
for  the  formula  of  R  in  terms  of  F).  The  following  routines  have  been 
developed  and  used  in  a  variety  of  problems. 

1 .  Multidimensional  least  square  technique. 

2.  "Overlapping"  least  square  technique. 

3.  Fourier  decomposition  of  each  section  of  a  surface  and  then 

blending  them  to  obtain  the  equation  of  a  surface. 

Beside  the  development  of  the  above  techniques,  we  have  also  developed 
the  following  two  separate  programs  in  the  solution  of  Eq.  (21). 

1.  For  the  acceleration  of  the  iterative  process,  a  "multigrid"  tech¬ 
nique  in  the  solution  of  Eq.  (21)  has  been  incorporated.  From  the 
test  cases  conducted  so  far  it  looks  that  the  computation  time  with 
multigrid  is  definitely  much  lower  than  with  ordinary  SOR.  At  the 
time  of  the  progress  report  in  Agusut  1987,  these  results  were  not 


definite.  The  application  of  the  multigrid  technique  in  the  solu¬ 
tion  of  the  surface  grid  equations  seems  to  be  the  first  attempt  of 
its  kind. 

2.  A  routine  for  the  calculation  of  the  optimum  acceleration  parameter 
has  been  developed  which  works  quite  well  with  the  SOR  technique. 

3.  CONCLUSIONS 

The  problem  of  numerical  coordinate  generation  in  arbitrary  surfaces 
has  been  addressed  by  first  developing  a  mathematical  model  as  formed  of  a 
set  of  elliptic  PDE's  and  then  solving  the  proposed  equations  numerically. 
The  proposed  mathematical  model  has  deliberately  been  made  to  depend  on  the 
formulae  of  Gauss  which  involve  the  partial  derivatives  of  the  Cartesian 
coordinates  with  respect  to  the  curvilinear  coordinates.  With  the  formulae 
of  Gauss  as  the  basis  for  the  proposed  elliptic  equations,  it  is  logical  to 
conclude  that  on  solution  the  generated  Cartesian  coordinates  will  be  second 
order  differentiable  with  respect  to  the  curvilinear  coordinates. 

The  proposed  equations  have  been  used  to  generate  the  Cartesian  coordi¬ 
nates  as  functions  of  the  curvilienar  coordinates  in  a  variety  of  surfaces. 
The  forcing  functions  in  the  PDE's,  which  distinguish  one  surface  from  the 
other,  depend  on  the  equation  of  the  surface  in  the  form  F(x,y,z)=0.  When  a 
surface  is  given  only  through  a  discrete  set  of  data  points,  the  method  of 
least  squares  and  a  Fourier  method  have  been  developed  to  fit  the  required 
function  F. 

A  computer  program  has  been  developed  for  obtaining  the  curvilinear 
coord inates/grids  in  either  a  simply-  or  doubly-connected  type  of  surface. 
For  complicated  shapes,  a  segmentation  approach  is  used  in  which  the  whole 
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ABSTRACT 

This  paper  is  devoted  to  an  analytical  comparison  of  the  various  elliptic  partial-dif¬ 
ferential-equation  (PDE)  models  which  are  in  current  use  for  grid  generation.  These 
comparisons,  particularly  between  the  equations  from  the  Laplace-Poisson  system  and 
the  equations  from  a  Gaussian  approach,  have  yielded  useful  expressions  connecting 
the  3D  Laplacians  and  the  surface  Beltramians.  This  effort  has  specifically  been 
successful  when  the  transverse  coordinate  leaving  the  surface  is  orthogonal  to  the 
surface.  Equations  which  are  derivable  from  Cartesian-type  Poisson  equations  and 
those  obtained  by  using  the  variational  principle  in  surface  coordinates  have  also  been 
considered. 


I.  INTRODUCTION 

The  problem  of  generating  spatial  coordinates  by  numerical  methods  is  of 
much  importance  in  many  branches  of  engineering  mechanics.  A  review  of 
various  methods  of  coordinate  generation  in  both  two-  and  three-dimensional 
Euclidean  spaces  is  available  in  Thompson,  Warsi,  and  Mastin  [1,  2). 

This  paper  is  exclusively  directed  to  a  collection  and  analytical  compari¬ 
sons  of  the  various  elliptic  partial-differential-equation  (PDE)  models  which 
are  curr'ntiy  in  use  for  numerical  coordinate  generation. 

The  theory  of  grid  generation  does  not  depend  on  any  set  of  so-called 
conservation  laws,  and  thus  a  variety  of  equations  and  methods  of  different 
characters  can  be  used  to  obtain  the  grids.  Any  consistent  method,  depending 
either  on  the  solution  of  PDEs  or  any  algebraic  method,  can  be  used  to  obtain 
intersecting  trajectories  in  either  2D  or  3D  Euclidean  space.1 


“Research  supported  by  the  Air  Force  Office  of  Scientific  Research  under  Crant  AFOSR-85-0143 
‘The  question  of  space  comes  into  the  picture  when  it  is  realized  that  all  the  metric  coefficients 
g,(  cannot  be  selected  arbitrarily.  In  fact,  these  metric  coefficients  must  be  so  selected  that  a  set 
of  second-order  PDEs  become  satisfied.  For  details  refer  to  Warsi  [3|. 
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In  light  of  the  above  discussion,  any  consistent  set  of  PDEs  is  sufficient  to 
form  a  mathematical  model.  In  bounded  domains  the  most  natural  choice  is 
-  that  of  a  mathematical  model  formed  of  elliptic  PDEs.  The  simplest  set  of 
equations,  which  also  yields  the  smoothest  grid,  consists  of  the  Laplace 
equations  of  the  curvilinear  coordinates  in  the  Cartesian  physical  space.  As  a 
next  logical  step,  a  set  of  Poisson  equations  can  be  selected  so  as  to  have  a 
degree  of  control  of  the  distribution  of  grid  lines.  On  inverting  these  equations 
a  set  of  quasilinear  PDEs  are  obtained  [see  Equations  (2),(4)J. 

After  the  development  of  the  problem  of  grid  generation  through  solving 
the  inverted  forms  of  the  Laplace  and  Poisson  equations  in  2D  Euclidean 
domains  by  Allen  [4],  Winslow  [5],  Chu  [6],  and  Thompson  et  al.  [7],  a  logical 
extension  was  to  use  the  same  equations  for  3D  domains,  as  has  been  done  by 
Thompson  and  Mastin  [8].  Parallel  to  the  above-noted  developments,  Warsi 
[3,  9,  10]  proposed  a  Gaussian  approach  which  basically  generates  surfaces 
and  thus  can  be  used  either  for  generating  curvilinear  coordinates  in  a  given 
surface  [11]  or  for  generating  3D  coordinates  by  generating  a  series  of 
surfaces  starting  from  the  data  on  the  given  surfaces  [9,  12].  The  Gaussian 
approach  in  fact  depends  on  a  manipulation  of  the  formulae  of  Gauss  for  a 
surface,  and  thus  the  resulting  equations  have  the  surface  coordinates  as  the 
independent  variables.  This  manipulation  introduces  the  Beltramians  of  the 
curvilinear  coordinates  and  the  sum  of  the  principal  curvatures  of  the  surface 
in  a  natural  way.  Since  the  formulae  of  Gauss  for  a  surface  hold  true  for  any 
allowable  coordinate  system  introduced  in  the  surface,  the  equations  pro¬ 
posed  by  Warsi  [3,  9,  10]  must  also  have  the  same  properties.  Further, 
because  of  the  use  of  Gauss  formulae,  the  proposed  equations  are  optimal  in 
the  sense  that  the  number  of  terms  and  the  amount  of  information  in  the 
equations  is  just  sufficient  for  the  generation  of  either  surfaces  or  coordinates. 

In  this  regard  we  can  justifiably  call  a  mathematical  model  '‘optimal’'  if  it 
can  be  reduced  to  the  form  of  the  proposed  equations  [cf.  Equation  (19a)]. 

Beside  the  Laplace-Poisson  system  and  the  Gaussian  approach,  we  have 
also  derived  the  surface  generating  equations  by  the  use  of  the  variational 
principle.  The  resulting  equations  are  near  optimal  for  the  Gaussian  system. 


2.  BASIC  ELLIPTIC  MODELS 

Poisson  Equations  as  Grid  Generators 

Since  the  publication  of  the  TTM  method  [7],  there  has  been  extensive  use 
of  the  Poisson  equations  in  the  physical  r-space  to  generate  both  2D  and  3D 
grids  [2].  In  practically  all  cases  the  main  aim  is  to  have  those  equations  in 
which  the  computational  coordinates  appear  as  the  independent  variables. 
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and  therefore,  the  Poisson  equations  have  to  be  inverted  by  making  the 
physical  coordinates  r  =  (x,  y,  z)  the  dependent  variables.  This  inversion  can 
of  course  be  carried  out  in  a  nontensorial  manner  by  using  the  chain  rule  of 
partial  differentiation.  However,  it  is  much  simpler  to  follow  either  of  the 
following  two  methods  to  attain  the  same  result. 

Method  I:  Inner  multiplication  of  Equation  (A2)2  by  g‘>  and  the  use  of 
Equation  (A4)  results  in 


(1) 


Introducing  the  second-order  differential  operator 


D  =  ggij 


d2 

ax’ dx> 


into  Equation  (I),  we  get 


Dr+  g(v2xk)r  k  =  0,  (2) 

which  is  the  desired  equation  in  vector  form  in  the  computational  plane. 

Method  II:  In  Equation  (A3),  writing  <f>  =  r  =(x,y,z)  and  again  using 
Equation  (A4),  we  obtain  the  same  equation  (2). 

To  form  a  closed  system  of  equations  from  Equation  (2)  one  has  to  specify 
the  Laplacians  V2!*  in  a  suitable  fashion.  As  described  in  Reference  [3],  a 
general  specification  for  the  Laplacians  is  to  write 


v  V^g'jpk,  fc  =  1,2  or  k  =  1,2,3,  (3) 


and  then  Equation  (1)  takes  the  form 

g"(r,,,  +  ^r.ik)  =  0.  (4) 

In  Equation  (3),  the  P*  =  P-  are  intended  to  be  arbitrarily  specified  control 
functions.  Equating  the  right-hand  sides  of  Equation  (A4)  and  Equation  (3), 


2  Refer  to  the  Appendix. 


vv.v, . 


msmmmm 
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we  have 


g.ip*= -g.yr*.  (5) 

Since  g‘>  are  not  arbitrary,  we  conclude  that 

Pk*  -  r* 

>/  u 

for  all  values  of  i,  j,  1:. 

Recursive  Property  of  P* 

We  now  impose  the  following  requirement  on  Pk:  //  the  coordinates  are 
such  that  their  Laplacians  vanish,  i.e.,  Vzxk  —  0,  then  the  control  functions 
Pf  vanish  individually  for  all  indices  i ,  j,  k. 

The  importance  of  this  requirement  becomes  clear  by  considering  Equa¬ 
tion  (5)  in  which  T,k  are  the  Christoffel  symbols.  Thus,  when  the  Laplacians 
of  the  coordinates  vanish,  the  right-hand  side  of  Equation  (5/  vanishes  as  an 
inner  sum,  but  the  left-hand  side  vanishes  due  to  the  imposed  restriction  of 
“ individual  vanishing”  of  Pk  for  such  coordinates.  Denoting  the  coordinates 
which  satisfy  the  Laplace  equations  by  xj0),  we  then  have 

^<0,  =  °-  r*(0)*o,  g<V?;(0)  =  o.  (e) 

We  now  consider  those  coordinates  which  satisfy  the  Poisson  equations.  In 
these  coordinates  it  is  of  interest  to  know  the  relation  between  the  successive 
P  's  under  successive  coordinate  transformations.  To  this  end,  we  consider  two 
successive  allowable  transformations  denoted  by  x(‘m_  n  and  xjm),  with  xj0)  as 
those  coordinates  which  satisfy  the  Laplace  equation.  Thus,  in  the  x|m_  ^ 
coordinates.  Equation  (4)  is 


a  1 1 
o(  m  -  1) 


dzx 


and  in  the  rJTO)  coordinates 

d2r 


+  P.p. 


dr 


jj(m-  1>  t  o 


=  0, 


+  P*. 


=  0. 


(8) 
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Using  the  transformation  law  (AT)  in  the  form 


3x{m-  1 


k(m) 


(m ) 


2pn 
&(  m) » 


0) 


and  the  chain  rule  of  partial  differentiation  of  r  in  Equation  (7),  and 
comparing  the  resulting  equation  with  Equation  (8),  we  get 


<5xfm_n 


(m  -  l) 


3x(m-D  9x"m-l) 


dx! 


(m  ) 


3x'm) 


.  (10) 


establishing  a  relation  between  the  successive  P ’s.  This  is  how  P*  transform 
from  one  coordinate  system  to  the  other.  .Also,  using  Equation  (A9)  in 
Equation  (10),  we  have 


pp 
i  j(  m ) 


l'q'(m)  ^  l'lcn(m-ll) 


3x{m) 


dxtm-U  3x"m-l) 


dx 


lm-1) 


dx 


(m) 


3x[m) 


(ID 


Equation  (11)  establishes  a  relation  between  the  P’s  and  T’s.  In  particular, 
for  m  =  1,  Equation  (10)  reduces  to  Equation  (63)  of  Reference  [3). 

From  the  recursive  relation  (10)  it  is  a  straightforward  matter  to  show  that 
the  transformation  x(‘0)  -*  x^  x(2)  of  P,k  is  the  same  as  the  transformation 
x(0)  x(2)'  The  ^  °1  tl|e  chain  rule  for 


3x(2)  32x’(2) 

3xtl)  ’  dx(l)  3xiU 

and  of  the  formula 

fcfo)  3xw  dx\o)  dxm 

9x(l)  9xw  9xltn  9xw  dxii)  dx'u  9xii) 
gives  the  required  result. 

Other  Poisson  Systems 

It  is  also  possible  to  have  equations  in  which  the  dependent  variables  are 
non-Cartesian,  e.g.,  cylindrical  or  spherical  coordinates.  In  a  paper  by  Ghia 
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et  al.  [13],  the  computational  coordinates  are  assumed  to  satisfy  a  Poisson 
system  in  the  cylindrical  coordinates. 

To  generalize  this  concept,  let  x‘  be  a  coordinate  system  on  which  the 
coordinate  system  x1  (e.g.,  cylindrical)  is  to  be  generated.  From  Equation 
(A4), 


VI 


’  =  -  grip*. 


(12) 


and 


VV=  -g"T*  =  g  "P* 

=  F*  (say).  (13) 

In  Equation  (12)  both  the  metric  coefficients  and  the  Christoff  el  symbols  are 
already  known  in  terms  of  the  x ’  coordinate  system.  [If  x’  are  the  Cartesian 
coordinates,  then  Equation  (12)  is  an  identity.]  From  Equation  (A3)  we  then 
have 


g 


a *t* 
dx'dxi 


(14) 


which  are  the  transformed  equations  in  the  computational  space.  By  using 
Equation  (A10)  we  can  also  write  Equation  (14)  as 


gmnc;q 


a2?1  dx’  2- 

7  +  /  p  I?"-(/)  «',r‘V 


n  dx'dx 


(15) 


where 


dx* 


Recendy  Fujii  [14]  used  a  Cartesian-type  Poisson  system  between  two 
curvilinear  coordinate  systems.  Let  x5  and  x’  be  two  curvilinear  systems. 
Consider  the  Cartesian-type  Poisson  system 


I 


aV 

dxpdxp 


=  p\ 


(16) 
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where  x'  is  the  curvilinear  system  on  which  the  xp -system  (e.g.,  spherical)  is 
.to  be  generated.  To  obtain  the  inversion  of  Equation  (16),  we  take  the  inner 
product  of  Equation  (All)  with  dx'/dx'  to  obtain 


d2x‘  dx'  3x> 
dx‘  dx’J^lP 


<?V  dx1 


3xkdin  dxr' 

Setting  k  =  n  =  1,  k  =  n  =  2,  k  -  n  =  3  and  adding,  we  get 


C> 


d%*  dx'  _  dV 


dx‘dx> 


>  dx,/Ljaz 


dxp  dxp' 


(i‘ 


so  that 


C‘i 


d2x‘  dx' 


dx'  dx 


-  = - P' 

i  dxr  ’ 


where 


(18) 


dx 1  dx>  1  _ 

Equation  (18)  is  the  inverted  form  of  Equation  (16).  Similarly  to  the  case  of 
actual  Poisson  system,  the  functions  PT  are  again  the  arbitrary  specified 
functions.  Alternatively,  we  can  also  write 

P'  =  C''P'j, 

and  then  P'-  are  the  arbitrary  control  functions. 

Beltramians  as  Grid  Generators  in  Curved  Surfaces 

The  problem  of  generating  curvilinear  coordinates  in  a  surface  has  im¬ 
portant  applications  in  many  branches  of  engineering.  In  this  paper  we  are 
interested  only  in  one  differential  model  which  generates  coordinates  in 
surfaces  embedded  in  R  3  and  which  reduces  to  Equation  (2)  in  a  natural  way 
when  the  surface  degenerates  into  a  plane.  To  achieve  this  aim,  we  consider 
the  formulae  of  Gauss  [Equation  (A12)]  for  a  surface  v  =  const.  Inner 
multiplication  of  Equation  (A12)  by  C„ gaS  yields 


Dr+G,(A(;'is)r4  =  n(')R, 


(19a) 


rj 

tj 

s 
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3.  SURFACE-GRID  EQUATIONS  FROM  THE  LAPLACIANS 

In  this  section  an  attempt  is  made  to  establish  a  relation  between  the 
surface  equations  derived  by  using  the  formulae  of  Gauss  [Equation  (19a)] 
and  the  surface  equations  as  derived  from  the  inversions  of  the  3D  Laplacians. 

Let  t),  f  be  a  general  curvilinear  coordinate  system  in  R3  such  that  £,  7) 
form  the  coordinates  in  the  surface  f  =  const,  with  f  as  the  transverse 
coordinate.  Starting  from  Equation  (B.3)  of  Reference  [1]  and  using  Equation 
(A2),  we  get  (refer  to  [11]  for  more  details) 

+  Mi1  )rf 


where  the  differential  operator  (2  is 


fl  =  D  +  T° 


~D  +  (S“  +  C3  Asx“)- 


In  (22b,  c)  the  differential  operator  D  defined  in  (19b)  is 
D  ~  g»2^££  —  "**  Sll^rn’ 


2'“  — 2gl2ri‘,2  gzoTj”  gu^m> 

sa  =  2g„(rs  -  Tjk)  -  gjrz  -  TfJ  -  gu(r&  -  t&). 


where  a  =  1,2.  Further 


X  =  n-rr. 


n  being  the  unit  surface  normal  vector  on  f  =  const. 

It  must  be  emphasized  that  Th  and  are  the  3-space  and  2-space 
Christoffel  symbols  of  the  second  Kind  respectively,  and  in  general  they  are 
not  equal  to  each  other  at  f  =  const.  Thus,  a  comparison  of  Equations  (19a) 
and  (22a)  shows  that  the  two  sets  of  surface  generators  are  entirely  different. 
However,  it  has  been  shown  below  that  Equations  (19a)  and  (22a)  become 
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exactly  the  same  equations  under  the  following  conditions: 

(i)  When  the  surface  degenerates  into  a  plane,  in  which  case  £,  rj  are  the 
curvilinear  coordinates  (of  any  nature)  in  a  plane. 

(ii)  When  £,  -q  are  ?ny  general  coordinates  in  a  surface,  but  f  (the 
transverse  coordinate)  is  orthogonal  to  the  surface. 

Case  (i)  is  patently  straightforward.  For  case  (ii),  noting  that  f  is  orthogo¬ 
nal  to  the  surface,  we  have 


8l3  — 823  0* 

(23a) 

and 

n—  I -  ^  —  1/833  • 

V633  > 

(23b) 

Under  the  conditions  (23a)  it  is  easy  to  show  that 

11 

£ 

where  all  Greek  indices  assume  only  two  values  (here  1,2).  Thus  S“  =  0,  and 
using  (23b)  we  find  that  Equations  (19a)  and  (23a)  become  the  same 
equations  under  the  condition  of  orthogonality  of  the  f -coordinate  to  the 
surface  f  =  const.  Further,  under  this  condition  the  Beltramians  and  the 
Laplacians  of  $  and  tj  are  related  through  the  following  equations: 

r1 

2>  a  j.  A33 

V2$  =  A2£ - , 

833 

(24a) 

r2 

2  a  *33 

V  7j  =  A  ,77 - - 

833 

(24b) 

Further, 

r2,  k?  +  k<? 
v  -  —  , 

633  V833 

(24c) 

where 
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and  [referring  to  (19b)] 

^3  =  SnSis  ~  (&12)  • 


4  SURFACE  GENERATING  EQUATIONS  BASED  ON  A 
VARIATIONAL  PRINCIPLE 

The  use  of  a  variational  principle  in  the  generation  of  both  2D  and  3D 
grids  has  already  been  considered  by  Brackbill  and  Saltzman  [15,  16]  and 
Thompson  and  Warsi  [17].  In  this  section  we  shall  consider  only  the  surface 
grid  generation  problem  based  on  the  use  of  a  variational  principle.  In  essence 
the  following  analysis  is  a  unified  approach  to  both  the  plane  2D  and  surface 
2D  cases  from  a  variational  viewpoint. 

Let  x1,  x2  be  the  coordinates  in  a  surface.  Consider  the  surface  functional 

Z  =  fjc^tdx'dx2,  (25) 

where  G3  =  gngn-(gl2)2  and  <f>  is  a  specified  function.  The  condition 
SI  =  0  then  leads  one  to  the  Euler  Lagrange  equations  (using  the  summation 
convention) 


3  sjjclfl 

dxff  dxr  p 


=  0, 


(26) 


where  xr  (r  =  1,2,3)  are  the  rectangular  Cartesian  coordinates,  x^  are  the 
curvilinear  surface  coordinates,  and 


dxr  =  d2!, 

dx“’  dxadx^~ 

From  Equation  (26),  it  is  a  purely  algebraic  problem  to  show  that 

JG~2  d  (  l  <?C3  \ 


(27a) 


where 


Let  now  <p  be  a  function  of  xr  p.  Then  expanding  Equation  (26)  and  using 
Equation  (27),  we  get 


Lx.  =  -  M, 


where 


ac3 

dip  dip 

dC3 

3 

d 

dip  \ 

dxr.fi 

dxfi  +  dx, p 

dx\ 

<p 

dx 2 

dXr.0l 

<t> 

M  =  - 
2 


On  the  other  hand,  if  <j>  is  taken  as 

<p  =  f/c3, 

where  F  is  still  a  function  of  xr  p,  then  in  place  of  (28b),  we  have 

1 


(28a) 


(28b) 


(29a) 


M  = 


2  F 


dC3  8F  dC3  dF 


dxr  p  dxp  dx^  dxT  B 


1  dC3  dC3  G3  d 
2GZ  dxr,plP~TlP 


3F 

\  dXT.0 


(29b) 


The  generating  system  (28a)  with  L  defined  in  (27b)  is  similar  to  Equation 
(19a).  However,  the  selection  of  the  form  of  the  function  <j>  or  F  which  yields 
the  right-hand  side  of  Equation  (19a)  seems  to  be  a  difficult  task.  One  simple 
case  is  when  tp  =  1.  In  this  case  the  minimization  of  I  implies  [from  (28a)] 

Lx,  =  0, 

and  these  are  the  equations  for  a  minimal  surface.  The  form  (29a)  is  of 
interest  because  the  choice 


F  —  f>u"^S'22>  (oOa) 

or 

<*>  =  gH  +  g-,  (30b) 


is  equivalent  to  the  "smoothness"  problem  in  2D  plane  coordinates,  as  has 


i  m  n  iiH  u  m  wm  ii»  ■  *irvir\ 
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been  shown  in  Reference  [15].  It  must,  however,  be  stated  that  "smoothness” 
of  coordinates  in  a  2D  plane  problem  is  due  to  the  satisfaction  of  the  Laplace 
equation.  No  such  criterion  is  obvious  on  using  either  of  the  equations  (30)  in 
(28a). 


5.  ANALYTICAL  COMPARISONS  AND  CONCLUSIONS 

Based  on  the  foregoing  analysis,  we  conclude  as  follows: 

(a)  The  Laplace-Poisson  system  for  2D  regions  is  optimal,  since  its  inver¬ 
sion  coincides  exactiy  with  the  Gaussian  equations  in  a  plane.  [Note:  In  a 
plane  the  right-hand  side  of  Equation  (19a)  is  zero  and  the  Beltramians 
become  the  Laplacians.] 

(b)  The  inversion  of  the  Laplace-Poisson  system  for  3D  regions  for  a 
constant  coordinate  value,  viz.,  for  a  surface,  does  not  reduce  to  the  Gaussian 
equations  (cf.  Section  3)  except  when  the  transverse  coordinate  is  taken  as 
orthogonal  to  the  surface  [cf.  Equation  (22a)].  This  implies  that  the  Laplace- 
Poisson  system  in  3D  regions  with  three  nonorthogonal  coordinates  is  not 
optimal,  though  it  is  a  valid  system.  The  extra  terms  (22f)  should  somehow  be 
managed,  and  in  practice,  they  are  taken  as  part  of  the  arbitrary  specified 
control  functions.  This  means  that  the  generated  coordinates  will  assume  a 
distribution  which  may  not  be  to  one’s  desire.  It  must,  however,  be  again 
emphasized  that  the  terms  (22f)  vanish  when  the  transverse  coordinate  is 
orthogonal  to  the  surface.  (In  Reference  [18],  the  author  had  to  make  other 
assumptions  besides  orthogonality.) 

(c)  For  the  elliptic  system  described  by  Equation  (15)  the  conclusions 
discussed  in  (b)  hold  good. 

(d)  The  generating  system  described  by  Equation  (16)  is  a  Cartesian-type 
Poisson  system.  It  looks  difficult  to  assess  its  optimality  in  relation  to  the 
Gaussian  equations. 

(e)  The  equations  derived  from  the  variational  principle,  viz.,  (28),  are 
nearly  optimal,  though  it  looks  difficult  to  find  the  appropriate  function  <t> 
which  makes  the  right  hand  side  of  (28a)  the  same  as  that  of  Equation  (19a). 


APPENDIX 

The  following  formulae  have  been  used  in  the  main  text  of  the  paper  and 
can  be  found  in  any  standard  text  on  tensors,  e.g.,  [19],  [20]. 

In  an  Euclidean  space  R "  (though  here  we  are  concerned  only  with  R -  or 
R3),  endowed  with  the  curvilinear  coordinates  x',  the  first  partial  derivatives 
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of  the  covariant  basis  vectors  a,  are  expressible  as  linear  functions  of  a,: 


fa. 

dx ’ 


r‘a 


ijuk* 


(Al) 


where  T,*  are  the  Christoffel  symbols  of  the  second  kind.  Here,  and  in  what 
follows,  a  repeated  lower  and  upper  index  will  always  imply  summation  over 
the  range  of  index  values.  Equation  (Al)  can  alternatively  be  written  as 


r.if  = 


(A 2) 


where  a  comma  preceding  an  index  implies  a  partial  derivative,  and  r  is  the 
position  vector,  viz.. 


r  =  (x,t/,=)or(x1,i2,x3). 

The  Laplacian  of  a  scalar  <p  is  given  by 

V24>  =  gi,UJ/-r*4>,k).  (A3) 

From  (A3),  we  have 

V  V  -  -  giJT,i'  (A4) 

The  quantities  g‘J  and  g,;  are  respectively  the  contravariant  and  co variant 
components  of  the  metric  tensor,  and  the  two  are  related  as 


where  (t,  p,r)  and  (;',  n,s)  are  to  be  taken  in  the  cyclic  permutations  of 
(1,2,3),  and 


g  =  det(g, .). 

Let  x‘  and  f  be  two  allowable  coordinate  systems  in  a  Euclidean  space 
such  that  each  of  the  functions  x'  =  <J>'(x J)  and  x‘  =  define  a  transfor¬ 

mation  with 
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On  transformation  of  coordinates  from  r‘  to  x',  we  have 


g  j  2  •  g  j2  ’ 

dx'  dx> 

gpn  ~ 

dxp  di" 

6  dx'  dxi°  ' 

_  dxp  dx'  dx ’  d^x’  dx” 

kn  '<  dx’  dxk  dxn  dxkdxn  dx>’ 

d2xp  ^  dxp  _n  dxr  dx' 
dxkd~xn  =  Vkn  a?  “  r,‘  a?  aT" ' 

The  first  partial  derivatives  of  r'  with  respect  to  f  ’  are  given  by 

dr*  C! 


dx'  J  ’ 


where 


dxr  dxk  dxT  dxk 

(AI0) 

and  (i,  s,  n )  and  (;',  r,  k)  are  cyclic  permutations  of  (1,2,3V 

The  second  partial  derivatives  of  one  set  of  coordinates  are  related  with 
those  of  the  other  set  as 


d2v  dxr  dx ■  dx' 


dik  dx” 


dx'dx >  dx’  dxk  d£n ' 


We  now  consider  a  surface  embedded  in  R3.  All  tensor  indices  associated 
with  a  surface  will  be  denoted  by  Greek  letters  (except  the  letter  r).  In 
contrast  to  Equation  (A2),  the  formula  of  Gauss  is 

r.afl  =  TaV,5+n(',ka<i.  (M2) 
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where  are  the  'urface  Christoffei  symbols  of  the  second  kind,  baB  are  the 
coefficients  of  the  second  fundamental  form,  and  n(,)  is  the  unit  surface 
normal  on  the  surface  v  =  const.  The  values  of  v  and  other  Creek  indices 
follow  the  following  scheme: 

v  =  1:  Greek  indices  a,  0,  etc.  assume  values  2  and  3. 

v  =  2:  Creek  indices  a,  0,  etc.  assume  values  3  and  1. 

v  =  3:  Creek  indices  a,  0,  etc.  assume  values  1  and  2. 

The  Beltramian  of  a  scalar  <p  is  given  by 

A(r*  =  gaVa*-T>.5).  (A13) 

From  (A13),  we  have 

A';'xs=  -ga07sa$.  (A14) 

For  a  surface  the  formulae  (A5)-(A9)  are  equally  applicable  with  proper 
choice  of  indices  and  replacing  T  by  T. 
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